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Abstract 

In this work we derive the evolution equation for the den- 
sity contrast considering open system cosmology, where the 
influence of adiabatic particle production process on the dy- 
namic of a homogeneous and isotropic is investigated within a 
manifestly covariant formulation. As application we derive the 
solution for two sources, one of them is a generalization of Pri- 
gogine's model. Then we establish the condition for the reach 
of the non-linear regime for the density contrast which turns 
out to be a necessary condition for the structure formation. 



1 Introduction 



The mathematical theory of perturbations in homogeneous, isotropic 
models has been worked over many times in the literature [[TJ]- 0. 
Nevertheless, troubling questions still remain about the physical 
interpretation of density perturbations at early times when the per- 
turbation is larger than the particle horizon, which will here mean 
when the time for light to travel a characteristic wavelength of 
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the perturbation is larger than the instantaneous expansion time. 
These questions are particularly relevant to attempt to explain the 
origin of perturbations which eventually give rise to galaxy forma- 
tion. 

The study of perturbations of an expanding universe would ap- 
pear to be hopelessly gauge dependent. However, Bardeen [g] gave 
an important contribution to define perturbations that are gauge 
invariant non-geometrical quantities denned with respect to a par- 
ticular chart. Thus a gauge-independent theory of cosmological 
perturbations is possible. However, Bardeen's paper is difficult to 
understand and to apply. 

Ellis and Bruni gave an alternative approach to gauge invari- 
ant cosmological perturbations with basic variables more closely 
related to the physical quantities. They look for variables which 
vanish in the background. Quantities of this kind are always gauge 
independents. The gauge invariant key variable is the comoving 
fractional spatial gradient of the energy density. 

The purpose of this paper is to discuss an approach to cosmo- 
logical perturbations which are thermodynamically oriented gauge- 
invariant description of scalar cosmological perturbations around a 
flat homogeneous universe with particle production. All basic vari- 
ables have their physical meaning on "comoving hyper-surfaces" 
, orthogonal to the matter world lines. Thermodynamic perturba- 
tions relations in Minkowski space-time turn-out to be valid in the 
expanding universe, provided all perturbed quantities are replaced 
by their gauge-invariant counterparts. 

The characteristic feature of this model is the use of an explicit 
phenomenological balance law for the particle number with a non- 
zero source term together with the equation for the cosmic scale 
factor in a perfect fluid cosmology that includes a pressure creation 
that behaves like a viscous pressure, responsible for the particle 
production. 

2 Open system thermodynamics 

In the presence of matter creation, the appropriate analysis is per- 
formed in the context of thermodynamics of open systems [g] . This 
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consideration leads to an extension of thermodynamics as associ- 
ated with cosmology to reinterpretation in Einstein's equations of 
the matter stress energy tensor. In the case of adiabatic transfor- 
mation, the thermodynamical energy conservation leads to 

d{pV ) + P tk dV-^l = 0, (1) 

n 

where n = y is the particle number density and E = p + P th is the 
enthalpy per unit volume. N is the particle number and P th is the 
thermodynamical pressure. In such a transformation, the "heat" 
received by the system is due entirely to the change of the number 
of particles. In our cosmological context, this change is due to the 
transfer of energy from gravitational field. Hence, the creation of 
matter acts as a source of internal energy. 

The cosmic medium will be described by the energy momentum 
tensor 

T al3 = e^u 13 + Ph al3 , (2) 

where e is the energy-density, h a ' 3 is the projection tensor and u a is 
the four -velocity. P is the total pressure that includes the thermo- 
dynamical pressure and the creation pressure 0, namely, 

P = P th + P. (3) 

The creation process will be considered adiabatic, which here means 
a constant entropy per particle. Under this condition the equilib- 
rium entropy per particle does not change as it does in dissipative 
process. Instead, one can associate a "viscous pressure" to the 
particle production source, namely |9p 



n9 

where \& is the rate of particle production. 

Einstein's field equations imply the equations of motion 



and the energy balance 



h^ = 0, (5) 
u p T.f = . (6) 
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Additionally, according to the second law of thermodynamics, 
the only particle number variations admitted are such that 



h + n9 = m . (7) 

The general problem is split into a zeroth-order background uni- 
verse and first order perturbation about this background. The 
zeroth-order is assumed to be homogeneous, isotropic and spa- 
tially flat. The corresponding line element is 

ds 2 = -dt 2 + R{t) 2 ~i ab dx a dx b , (8) 

with the cosmic scale factor R(t) as the only relevant metric quan- 
tity. Consequently, Einstein's field equations reduce to 

3(|) 2 = e (9) 

2(f)- = -(e + P). (10) 

The dot denotes the derivative with respect to time. The surfaces 
t = const, are the space like surfaces of constant curvature and ^ ab 
is the 3-space metric. Consequently, the energy balance equation 
(6) becomes 

e = -3^(e + P). (11) 
Jri 

Considering the particle number density and the temperature as 
are our basic variables we can infer the state equations 

P th = P th (n,T) (12) 
e = e(n,T) (13) 
a = a(n,T), (14) 

where a is the entropy per particle. Making use of Gibbs relation 
and equation (14), it is straightforward to find the general relation 

de_ e + P th TdP th 
dn n n dT 
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Equation for evolution of the temperature is calculated considering 
relations (11) and (15), namely 

We can define the expression for the sound velocity as: 

= ^ (17) 
that can be reformulated using equation (15): 

\9P th l 2 
I dT J 

e + P th dn ' (e + P th )§ 
Taking into account the relation 



7 , 2 _ n dP th r[Vj nm 



Pth = ^ n+ ~dT T U9) 
and equations (7), (16) and (17), evolution of the pressure reads 

PtH =-^ + PtH){v 2 s + ^ r ——}-{e + P th ) (20) 
R ^ e + P th n 

dPth 

{ vl-^-}. (21) 

dT 

Finally the time evolution of the entropy density can be found by 
using Gibbs relation and equations (7) and (16), namely 

nTa = € -±p^{3n^z - ^} (22) 

where 

P 

z = — . (23) 

e + Pt h 

This relation provides a direct connection between the particle pro- 
duction rate and the expansion of the Universe. Considering the 
process as adiabatic together with (7) yields 

N 

z=f. (24) 



We sum up this section by writing the time evolution of the particle 
density, temperature, energy density and pressure in the adiabatic 
case, it is done by using (7), (16) and (19) 

n=-3|(l-z), (25) 



dPth 



3^T(^L) , (26) 



<9T 



e=-3^(e + P), (27) 



^th = -r-(e + P)v 2 s . (28) 



3 The perturbative equations 

We write the metric as 

ds 2 = (g a p + 5g afS )dx a dxP , (29) 

where g a/3 is the background metric tensor and 

5g ab = A5 ab + B Ab , (30) 
8gaA = F A , (31) 
5gu = E. (32) 

The quantities A,B,F e E are space and time dependent scalars 
hence we are interested in scalar perturbations. 

From normalization condition, = 1 and the line element (29) 
the four velocity in first order is given by 

R 2 E 

u A = u 4 = , (33) 

R 2 u° + R 2 F a = u a = v a , (34) 

where the hat indicates a perturbative quantities. Perturbing to 
first order the equation which gives the conservation of particle 
density (nu a ). a = ^, the result is 

As „ ■>„ / * R 2 E X 
- ■ + i?- 2 V V -VF+^ = - + -— . (35) 

n 2 n n n 2 



The energy balance and momentum balance equations to first 
order are respectively: 

\ + 3^(e + P) + (R- 2 Vv + { - VF)(e + P) = , (36) 

+ Pv,„ + (e + P)^ - R 2 ^^E ifl = . (37) 

To complete the set of basic equations we need to linearize field 
equations, namely 0: 

G\ = R~ 2 VA -^f + 2^VF - 3R 2 E = -e , (38) 

it it 



G 4a = _i a _ RRE ia - (R 2 + 2RR)F a (39) 
= -(e + P)v >a + R 2 PF, a . 

4 Density contrast evolution equation 

Variables representing perturbations introduced in the preceding 
section modify their values under gauge transformations. In the 
linear perturbation theory it is necessary only to consider infinites- 
imal transformation, namely 

x 'a = x *_£* (40 ) 

Zimdahl shows in his work [^] how to describe cosmological pertur- 
bations thermodynamically oriented by suitable quantities which 
are invariant under transformations (40). These objects have an 
obvious physical meaning and obey reasonable equations. The be- 
havior of scalars under transformation (40) follow 

S'(x)-S(x)=CS, n . (41) 

Therefore the scalar quantities of the preceding section transform 
as 

A' = A + 2^° (42) 
R 



7 



B = B + 2R~ 2 i (43) 

F' = F + R- 2 Z° + (R- 2 0,o (44) 

E' = E- 2fT 2 £° (45) 

v = v - £° (46) 

ft = ft - 3^<°(1 - z) (47) 
R 



r=T-3~nrr(-^-)(l-z). (48) 



9T 

The basic variables that we shall take are 



n R 
and 

f R 



h R 

- - 3^(1 - z) (49) 



^ - 3~«-fHl - *) • (50) 

A gauge invariant gravitational field variable from the point of view 
of the matter [|6p is: 

q = l A + 3 Ji v - (51) 

It will turn out that the variables q , r and s are sufficient for 
the gauge -invariant description of the system of perturbation equa- 
tions of the preceding section. 

The energy density perturbation and pressure perturbation are 
gauge invariantly characterized by: 

e c = e-3^v(e + P) (52) 

P c = Pth-S^v(e + P)v 2 s . (53) 

We can re-write the equation (49) by using (50): 

e c ={e + P th )s + Ta, (54) 

where 

de dP th 

a = — s . (55) 

dT dT 
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Introducing the gauge invariant variables, the particle number 
balance (35) may be written as 

s +(1 — z)q + -7— — — ^^—q (56) 

1 ,q 2 f Re + P th 3RR 2q 

R 2 . . , 7, R R if 

To find the last result it is necessary to use the perturbative field 
equations (38), (39), the field equation 

(|)- = "^ 2 7(1 " z) (57) 

and assumes the spatial dependence e ikr with constant comoving 
wavector k for all first-order quantities. 

When equation (56) is used into the energy balance (36) and the 
pressure balance (37), they become respectively 

(aR 3 )'TR~ 3 + + t = _ 3Hz ^h _ ZHP_ 



e + Pth 3RR n e + P th e + P t 



th 



and 



sv 2 s -vz + 3Hz(l - z)v{\ + v 2 ) + P n + ^q(l-z) (59) 

e + Pth 3R 



s 

^ a 



§ e + Pt h 

Substituting equation (59) into (56) it becomes 

s -3Hsv 2 + 9H 2 z(l-z)v(^~v 2 )-^^- (60) 

2 e + P t h 



— a-^ — h -7/2 + 3Hzs = — 

e + P th % 2 e + P t/l n 



Equation (60) can be written in terms of gauge invariant variables 
using equation (55), namely: 

s -3H^—^-(P c -le c )-(l-z)^t (61) 
e + Pth 2 c) 1 >3RR 
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We next focus our attention on the adiabatic case that is equivalent 
to consider a = [|5p, consequently equation (61) assumes the form 



2 7^ H . A hfq 



s - 3H(1 - z)(vi - -)s - (1 - z) — -V + 3Hzs = (62) 
2 J J 3RR 



where the variable s is given by 



e + Pfh 

A straightforward manipulation of equation (62) gives 



(63) 



S H-\{m-z){2vl- 1 )} + 3Hz+{Z + ^) + -^-}s (64) 

3 H7 
+ {-UH{l-z){2vl- 1 )}+3{Hz)+{3, - 



2 L V S 

- ^(2^- 7 )}{(l- 2 )(| + |) + i} + ^}5 

k 2 

- T2 vz = 0. 

This equation governs the behavior of the gauge invariant variable 
s that is proportional to the density contrast, in an universe with 
particle production at the expenses of the gravitational field. 

We recover the well known time evolution for the density con- 
trast when the particle production is absent, namely 

s+ ^-^ s = - 1651 

2 

where the scale factor of the background, RocU were used. 



5 Cosmological environment for particle pro- 
duction 

To integrate equation (64) the scale factor R(t) must be specified. 
This is done when the source of the particle production \[> is chosen. 
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Next, we study the time evolution for the perturbation of the 
matter density taking the source of the particle production given 

2 

by ^ = 3nj3H and ^ = aH » . Our purpose is to obtain the conditions 
for the reach of the non-linear regime in a "viscous universe" with 
particle production. 



5. 1 Source proportional to the product of the Hub- 
ble function and the particle number density 

Lima et al. [ fTQ] ] gave a reasonable justification to pick up ^ as 
^ = 3n/3H , where (5 is a positive constant . We can provide explicit 
solution for the time dependence of the expansion factor in the 
presence of particle creation by using the continuity equation 



state equation 

and the field equation 

where 



n + 3Hn = * , (66) 
P th = {y- l)p (67) 

RR + A(k + R 2 ) = 0, (68) 



A = |i/(l-/3)-l. (69) 



We obtain, taking k = and (3^1: 

R = R (h^. (70) 
h 

The subscript alludes to present time. Note that for (3 = 1 the 
dominant energy condition is violated [[TTJ . 

Substituting the scale factor (70) and the respective source of 
particle production in equation (64) , we obtain 

s + ds-C s = 0, (71) 

where 

6P{1 — p)t 
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and 

—2 4 A3 

C ° = -l^-3^W {3/? - 1+4/3 -V>- (73) 

To solve equation (71) we neglect terms that contain the sound 
velocity and consider the limit of the large wavelengths. Solutions 
for equation (71) are 



s_{t) = (74) 
s+ (t) = C 2 t^^ +2 . (75) 

In the absence of particle production, (3 = 0, relations (74) and 
(75) become the usual density contrast for Friedmann models with- 
out creation [ |12[ J . Although the decreasing mode can be important 



in some circumstances, we shall hereafter mainly deal only with 
the increasing mode. It is responsible for the formation of cosmic 
structure in the gravitational stability picture. Taking into account 
the decreasing mode the universe would not have been homoge- 
neous in the past. Besides, a growing mode that starts to grow just 
after the end of radiation era, Peebles argued that the contribution 
of decaying solution must be negligible [ |13| J . 

It is possible for the large perturbations that we see today to 
have grown from the small perturbation present in decoupling time. 
Typically we have(^) dec < (1(T 2 - 1(T 3 ) JT§. 

We compare the reach of the non-linear regime for open system 
cosmology with the same quantity for standard model, using the 
quotient given by: 

e = ^ , (76) 

0+F 

where the coefficient F refers to the Friedmann model and c is rel- 
ative to the model with particle production. In usual Friedmann 
universe the reach of the non linear regime of the density contrast 
is realized in a relative recent time, comparing with the age for the 
universe. Consequently, a cosmological model with better charac- 
teristics to represent the universe must be satisfied if > 1 . Taking 
into account the growing mode (75) and the condition (76) they im- 
ply that (3 < 0, for v = 1. Then we conclude that the respective 
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cosmological model to the source \& = 3n(3H does not give a better 
results than the usual FRW model, since that for negative values 
for /?, the second law of thermodynamics is violated. 



5.2 Source proportional to the H» 

2 

The source ^ = aH~ is a natural generalization of the Prigogine's 
work [|| for any value of v, excepted zero, where v is a parameter 
of the state equation (67). 

The field equations for the homogeneous and isotropic space 
time when particle production is taking into account can be written 
as: 

P=\h\ (77) 

71 

p=-(p + P th ). (78) 
n 

Using the state equation (67), the above field equations and the 
considered source we find for the scale factor: 

R(t) = {1 + a{e bt - 1)}£ , (79) 

where 

3 v ~,kM x x 
a o 

b = ^(f)i- (811 

The initial particle number density is denoted by n and M is the 
mass of the particle produced. 

Substituting the scale factor in equation (64) and changing the 
variable of the coordinate system using the following transforma- 
tion 

£ = ae bt , (82) 



we find: 



es(t)" + ^){6 + 3u + Jf }s(t)' + {— (83) 

3z/ 4 + 1 — a 3(4 + 1 — a) 

4£ 3f a-1. (a-1) 2 , . . 

3 a-1 +4 + ^-4 +4 ^ ; } 8 t =0. 

v l v v 3z/(4 + 1 — aj^ 
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The derivatives are now taking in relation to the variable £. The 
solution of last equation is given by the hypergeometric functions: 

s(0 = {C^hl, h2, h3; -^t) + C 2 F 2 (M, h5, h6; -^—) (84) 



Va-l 

where: 



zt) m R 2 » (f + a-1) 



fti = —{3 i/ + V1080 v 2 + 540 i^ 3 + 528 v + 81 z/ 4 + 4 (85) 
6z/ 

+ 3 V9 z/ 4 + 56 v 2 + 36 z/ 3 + 4 + 8 z/ + Vl6 + 24 z/ + 81 z/ 2 } 



/i 2 = —{3 v - V1080 v 2 + 540 z/ 3 + 528 v + 81 z/ 4 + 4 (86) 
6f 

+ 3V9 z/ 4 + 56 z/ 2 + 36 z/ 3 + 4 + 8 v + Vl6 + 24 v + 81 z/ 2 } 

z/ + V9 z/ 4 + 56 v 2 + 36 z/ 3 + 4 + 8 v 
h = (87) 



h 4 = —{3 v - V1080 v 2 + 540 z/ 3 + 528 v + 81 z/ 4 + 4 (88) 
6z^ 

- 3 V9 z/ 4 + 56 z/ 2 + 36 z/ 3 + 4 + 8 z/ + Vl6 + 24 z/ + 81 z^ 2 } 

/i 5 = — {3z/ - V1080 v 2 + 540 z/ 3 + 528 v + 81 z/ 4 + 4 
6z^ 

- 3 V9 z/ 4 + 56 z/ 2 + 36 z/ 3 + 4 + 8 v (89) 
+ V16 + 24Z/ + 81Z/ 2 } 



V9 z/ 4 + 56 v 2 + 36 z/ 3 + 4 + 8 z/ 

h 6 — u . (90) 

z/ 

The coefficients of the hypergeometrics usually obey the follow- 
ing relations: 

h 4 = h 2 + l-h 3 , (91) 
hh = h + 1 - h 3 . (92) 
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These conditions are valid for any value of v. A third condition 



h = 2-h 6 . (93) 

is valid only for v = 1. By inspecting the argument of the hyper- 
geometric function, one is convinced very quickly that no easy an- 
alytical solution is possible. In order to get an insight about the 
behavior of the evolution of the density contrast, we consider the 
constant a « 1 , consequently: 

£_>£_l + a . (94) 

Then, the equation (83) assumes the form: 

es(0" + C(l + ^)(10 + 3i/)a(0' + 2^4^(0 = . (95) 

6P AV* 

Solution of the last equation is given by: 

s(f ) = di wl - w2 + C 2 C 1+w2 06) 

where 

9^ 2 + 10 + 30z/ 
wl = (97) 

bp 

and 

V81Z/ 4 + 1080z/ 2 + 540^ 3 + 4 + 528// 

w2 = (98) 

bp 

When t — > 0, then £ — > I. So using the initial condition s(0) = 0, 
the integration constants must obey the relation d = —C 2 . For 
mathematical convenience we take C\ = § . 

Taking the solution for s, considering the parameter t — bt and 
i/ = 1, the anisotropy of density contrast correspondent to the 
anisotropy of the cosmic background radiation of the microwaves 
[s < 10 3 ) [|14| will occur in the interval 



r d < 4.1 x 10" 5 , (99) 
for p = 1. When r assumes the value 

n « 7.3 x 10~ 2 , (100) 
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the density contrast reaches the non linear regime. Now if we take 
t instead of r and considering that the anisotropy of CMBR occurs 
at « 10 6 years, then the structure formation occurs at 

t < 1.8 x 10 9 ys. (101) 

The density contrast in the FRW standard model reach the non- 
linear regime at t < 3.16 x 10 9 ys, consequently in open system cos- 
mology, for the source considered, the structures at large scale 
appeared in the universe more soon than in the standard model. 

6 Conclusions 

We obtain the evolution equation for the density contrast with re- 
spective solutions in the open system cosmology (OSC). 

According to our results one may say that the OSC is an alter- 
native model to describe the universe. The inclusion of the creation 
parameters indicates an older universe [ [T5| ] and furnishes a nega- 
tive pressure that accelerate the expansion of the universe, accord- 
ing the recent results obtained by the observations of the super- 
nova. We analyze the conditions imposed on the creation parame- 
ters considering the reach of the non-linear regime of the growing 
mode in a suitable time. The reach of the non-linear regime in the 
standard model occurs in a more recent time than in the open sys- 
tem cosmology, considering the source of particle production given 
by lit 
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